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By means of neutron interferometry the s- 
wave neutron scattering length of the He 
nucleus was re-measured at the Institut 
Laue-Langevin (ILL). Using a skew sym- 
metrical perfect crystal Si-interferometer 
and a linear twin chamber cell, false phase 
shifts due to sample misalignment were 
reduced to a negligible level. Simulation 
calculations suggest an asymmetrically 
alternating measuring sequence in order to 
compensate for systematic errors caused 
by thermal phase drifts. There is evidence 
in the experiment's data that this proce- 
dure is indeed effective. The neutron 
refractive index in terms of Sears' exact 
expression for the scattering amplitude has 
been analyzed in order to evaluate the 
measured phase shifts. The resuh of our 
measurement, b'^ = (6.000 ± 0.009) fm. 



shows a deviation towards a greater value 
compared to the presently accepted value 
of ^c = (5.74 ± 0.07) fm, confirming the 
observation of the partner experiment at 
NIST. On the other hand, the results of 
both precision measurements at NIST and 
ILL exhibit a serious 12a (12 standard 
uncertainties) deviation, the reason for 
which is not clear yet. 

Key words: elastic neutron scattering; 
few body systems; neutron interferometry; 
properties of the "^He nucleus; real part of 
the coherent scattering length. 

Accepted: August 11, 2004 

Available online: http://www.nist.gov/jres 



1. Introduction 

The scattering length a is defined as the zero-energy 
limit of the scattering amplitude / Owing to the spin 
dependence of the underlying strong interaction, the 
neutron scattering length must in general also be spin 
dependent. The spin independent part of the quantity is 
referred to as the coherent scattering length a^. In order 
to expand the scattering formalism to absorption, the 
scattering length is made complex, a = a' - ia'\ For a 
homogeneous, isotropic, monatomic and dilute gas 
sample of particle number density p, neutron scattering 
can be described by the so-called coherent wave that 
scatters off a macroscopic optical potential. Since this 
effective potential is weak in general, one can define a 
neutron refractive index and treat the scattering prob- 
lem similarly to geometrical light optics. When 



describing the scattering off a nucleus of mass number 

A in the laboratory system, one often includes the kine- 

matical transformation from the centre of mass system 

A + \ 
by defining the bound scattering length b a. 

A 

Sears [1] has given an exact expression of the neutron 
refractive index in terms of the scattering amplitude. 
Taking into account that the scattering amplitude /and 
the scattering length a are proportional only in the zero- 
energy limit, but that in general higher-order terms in 
the neutron wave number k contribute, one arrives at 
the following expression for the real (nonabsorptive) 
part of the refractive index in terms of the bound coher- 
ent scattering length b^: 



\-^pbx\-2kb:+o{k')i 



(1) 
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Since scattering lengths are of the order of a few fem- 
tometers, the second term in the bracket is of the order 
of 1 0^ for thermal neutrons and the term in J^ can safe- 
ly be ignored for our purposes. 

In the following section, we describe an interfero- 
metric experiment that measures the phase shift 



A(p = {n-\)kd 



(2) 



that the partial wave function of a neutron of kinetic 
energy E^ = fi^l^llm experiences when passing through 
a sample of thickness d. 



2. Experiment 



On their optical path through the interferometer, the 
partial beams gather a phase difference A<p that causes 
modulations of the detectable intensity behind the 
instrument: 



^,H "^ 1 + ^^^ ^*P- 



(3) 



By rotating of a plane parallel Al-plate through the 
divergent partial neutron beams after the dividing plate 
of the interferometer, an additional relative phase Acp^i 
between the partial wave functions can be superim- 
posed. The total phase shift is then the sum of this 
phase shifter contribution, a sample contribution 
^<Psampie ^^^ ^^ additional, time-dependent instrumental 
intrinsic phase A<Pint: 



The experiment was performed in May 2002 at the 
S 1 8 CRG facility at ILL using the skew symmetrically 
shaped perfect crystal interferometer IFM4. The setup 
is essentially the same as described in [2], the most 
important changes being a massive water-cooled 
Helmholtz coil of 50 cm diameter with the interferom- 
eter at its centre and two layers of Mylar^ foil between 
the coils' body and the interferometer. The coil was ini- 
tially installed for polarized interferometry but turned 
out to stabilize the instrument's temperature consider- 
ably when operated. 



A<p = A(p^i+A(p^^^pi^+A(p.^,. 



(4) 



Modulating the total phase by means of the phase 
shifter enables one to extract both modulus and sign of 
^<Psampie + ^<Ant- The intrinsic phase shift of the instru- 
ment can be monitored when the phase shift through an 
identical but evacuated sample cell is measured period- 
ically. The skew symmetry of the instrument has the 
consequence that the partial beams are almost parallel 
in the volume where the sample is placed, see Fig. 1. 
Therefore, the aluminum windows of the gas container 





Fig. 1. The partial beams emerging from the interferometer's crystal plates (T, SI, S2, A) enclose the double Bragg angle 29^. In 
contrast, at the sample's location the beams propagate nearly parallel. In the so-called sample position (left), the sample of refractive 
index n is placed in one partial beam, while the other partial beam passes through the evacuated compensation chamber. The twin 
cell is then moved to the so-called reference position, where both partial beams pass through the empty chamber (right). Now the 
phase shifter (PS) is rotated one angular step and one measures first in the reference position and then in the sample position. After 
a further rotation of (PS) the succession starts over until sample and reference measurements for 32 phase shifter positions 8 have 
been performed. 



Certain commercial equipment, instruments, or materials are iden- 
tified in this paper to foster understanding. Such identification does 
not imply recommendation or endorsement by the National Institute 
of Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the pur- 
pose. 
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do not act as phase shifter and any misahgnment of the 
container is compensated. Further, misalignment of the 
cell to better than 2.5° (which is easily achieved using 
a theodolite) results in an increase of the effective sam- 
ple thickness of less than 10^. 

The temperature correlated phase drifts of the instru- 
ment, c.f Fig. 2, would lead to a systematic error in the 
phase shift if they were not compensated for. 
Simulation calculations suggest an asymmetrically 
alternating measuring sequence. By this we mean the 
following procedure: for a fixed phase shifter position, 
neutrons are counted in both cell positions, first in the 
sample and thereafter in the reference position, say. 
Then, without changing the cell position, the phase 
shifter is moved one step and held fixed again. The sub- 
sequent two measurements are then performed in oppo- 
site order, first when the cell is in the reference position 
and then in the sample position. After that, the phase 
shifter is rotated one step further and the measure- 
ments' order starts over until data for 32 phase shifter 
positions are recorded. Performing alternating meas- 
urements with and without sample in changing succes- 
sion for every phase shifter position guarantees that 
even phase drifts during some 10 s are compensated. 

Moreover, the simulation calculation indicates that 
harmonic fits of the phase have to be carried out rela- 
tive to the center of the data. In this case, the absolute 
difference between the real phase shift and the fit result 
can in principle be made arbitrarily small. 
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Fig. 2. The intrinsic pliase of the instrument is obviously tempera- 
ture correlated. Phase drift rates up to (j) = n/{l2 h) were observed. 



The data show deviations from the expected sine 
behavior, which cannot be accounted for by Poisson 
statistics, see Fig. 3a. Instead, the deviations are 
assumed to be the result of integrating over short-time 
phase fluctuations. In order to include this error contri- 
bution, counting residuals were transformed into phase 
residuals, and histogrammed as shown in Fig. 3b. 
Neither the distribution of phase residuals nor the his- 
togrammed data show any repetitive pattern, thereby 
confirming the statistical nature of the underlying 
mechanism. Moreover, the standard deviation and the 
Gaussian width of the frequency distributions of the 
residuals are consistent in all cases. Typical phase fluc- 
tuations are of the order of some few tenths of radians. 
When this "phase fluctuation" contribution is com- 
bined with counting statistics, the data is well described 
by cosines. The reduced chi-squared values near 1 do 
not indicate any difficulties with this error analysis pro- 
cedure. 

Since the relative intensities of the partial beams 
were not constant during the experiment, transmission 
data could not be used to extract the gas particle num- 
ber density from the known absorption cross section. 
Instead, the temperature and pressure of the sample 
were monitored using DKD [3] calibrated sensors, and 
the density was calculated from the state equation of 
real gases. 

The neutron wavelength was determined to be X = 
(1.910 ± 0.002) A by measuring the angular difference 
between the Bragg peaks of the dispersive and nondis- 
persive crystal orientations. Figure 4 shows the meas- 
ured phase shifts of the individual measurements as a 
function of particle number area density. In accordance 
with Eqs. (1) and (2), the data show a strictly linear 
behaviour. Dividing the slope of the linear regression 
by the neutron wavelength and correcting for terms in 
^cZ?^' yields our final resuh of Z?,:= (6.000 ± 0.009) fm. 

We stress that there were over 1 3 h of measurement 
interruption combined with cell removal, soldering and 
cell replacement between the data points marked A and 
B. While one normally would avoid such "heavy work" 
during interferometric experiments, at the level of our 
precision the data taken after point B do not show a sys- 
tematic shift compared to the data taken before point A. 
Thus, we believe that the indications of the simulation 
calculations are confirmed by experiment and that our 
phase drift compensation scheme indeed works very 
effectively. 
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Fig. 3. a) Raw data (bottom) and residuals of the individual data points (top). Since the residuals do not show any systemati- 
cal pattern, the width of their frequency distribution b) can be used as an estimate for counting uncertainty contributions through 
integrating over fast varying phase fluctuations. 



10-1 



-% 



< 

en 

I 

Oh 



2- 




Phase shift He-n, 2=1.910(2)A 



B 



V ■'■'■ 



' Linear regression: 
A(p=1.628(4) + /?^ L146(l) 10"^V 

z'- 1,001 



1 



6x10 



24 



Gas particles per unit area p ^[m ] 

Fig. 4. Observed phase shifts as a function of gas particle number 
density times sample thickness. The slope of the linear regression 
divided by the neutron wave length yields essentially the real part of 
coherent bound scattering length, h'^i - 2kb'^+ OQ?)]. 



using largely the same technique. Compared to the lat- 
est NIST result, there is a 12 a discrepancy in the 
results, which cannot be resolved. 
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3. Discussion 

Compared to the currently accepted value of b'^ = 
(5.74 =t 0,07) fm, that also was measured using neutron 
interferometry [4], our result shows a significant 4 a 
deviation towards a higher value. This observation is 
confirmed by an independent experiment at NIST [5] 
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